An analysis of the ^2 



V (or extremely anisotropic next-nearest-neighbor Heisenberg) 
model 
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The t2 — V model (involving next-nearest-neighbor hopping and nearest-neighbor repulsion) has 
been shown to depict the limiting case of strong electron-phonon coupling in a molecular chain with 
cooperative breathing mode [R. Pankaj and S. Yarlagadda, Phys. Rev. B 86, 035453 (2012)]. Our 
t2 — V model can be mapped onto an extremely anisotropic Heisenberg model (with next-nearest- 
neighbor XY interaction and nearest-neighbor ising interaction) and is of the form J±_ '^i{Sf_iS~'.'^-\' 
H.c.) -I- J|| 5^^ SiSi^i- Using finite size scaling, at non-half-fillings of the t2 — V model (or non-zero 
magnetizations of the spin model), we numerically obtain the critical repulsion for a quantum phase 
transition by using modified Lanczos method. During the transition, away from half-filling (zero 
magnetization), the system undergoes a striking discontinuous transition from a superfiuid to a 
supersolid [i.e., a superfiuid homogeneously coexisting with a period-doubling charge-density-wave 
(antiferromagnetic) state]. At half-filling, the charge-density- wave (Neel) state and the superfiuid 
state are mutually exclusive. We also derive microscopically, using Green's functions, the exact 
instability conditions in the two limiting cases of the t2 — V model for hard-core-bosons: the two- 
particle system and the half-filled system. We show explicitly that the critical repulsion Vc for the 
two-particle case is Vc/i2 = 4 for a ring of any size while for the half-filled case the Vc/t2 = 2\/2 in 
the thermodynamic limit. Our spin model correspondingly lends itself to exact instability solutions 
(by the Green's function method) in the two limiting cases of two-magnons (the non-trivial highest 
excited state) and Neel antiferromagnet ground state. 

PACS numbers: 71.10.Fd, 67.80.kb, 71.45.Lr, 71.38.-k 



I. INTRODUCTION 



The cooperation or competition between diag- 
onal charge/spin long range orders and super- 
fluid/supereonductor ofF-diagonal long range orders 
in electronic phases leading to coexistence or mutual 
exclusivity of these orders is a subject of active current 
interest. Coexistence of charge-density- wave (CDW) 
and superconductivity/superfluidity is manifested in 
various systems with differing dimensions: helium-4^, 
bismuthates (e.g., BaBiOs doped with K or P)^, non-iron 
based pnictides (e.g., SrPt2As2)^, layered dichalogenides 
(e.g., 2H-TaSe2, 2H-TaS2, and NbSea)^, quasi-one- 
dimensional trichalcogenide NbSegi^ and doped spin 
ladder cuprate Sri4Cu2404i— , etc. 



A number of materials indicate strong electron-phonon 
(e-ph) interactions along with the larger energy scale 
electron-electron (e-e) interactions. For instance, oxides 
such as cupratesSriS, manganiteaii"— , and bismuthates^ 
indicate strong e-ph coupling. The coaction of e-e and e- 
ph interactions in these correlated systems leads to coex- 
istence of or incompatibility between various phases such 
as superconductivity, CDW, etc. 

Many oxides have a perovskite structure with the for- 
mula ABO3; they have two adjacent BOq octahedra 
sharing an oxygen which leads to cooperative octahe- 
dral distortions. Examples of simple three-dimensional 
systems exhibiting cooperative breathing-mode (CBM) 
phenomena are the barium bismuthates (BaBiOs) where 
Bi — O bonds of adjacent octahedra differ by about 



10%ii. Next, copper oxide a^^'^^ can be modeled as a 
one-band two-dimensional CBM system with the copper 
onsite energy being modulated by the movement of ad- 
jacent oxygen closer to or further from il^i^ii^. Lastly, 
in manganite systems when C-type antiferromagnetism 
manifests [as in Lai-xSr^MnOs for 0.65 < a; < O.Qii], 
the (1^2 orbitals participate in the C-chain ordering; the 
ferromagnetic C-chain is a one-band (involving d^^ or- 
bitals) one-dimensional CBM system. 

Earlier, we showed that cooperative effects in a one- 
dimensional CBM system change its dominant transport 
from a nearest-neighbor (NN) hopping mechanism to 
that of a next-nearest-neighbor (NNN) hopping^. Ad- 
ditional NN particle repulsion (due to incompatibility of 
NN breathing- mode distortions) leads to the t2 — V model 
as the effective model for the ID CBM model. 

In this paper, upon tuning repulsion, we show that 
the t2 ^ V model undergoes a dramatic discontinuous 
transition from a superfiuid state to a supersolid state 
(involving superfluidity in one sub-lattice only) . We find 
that (instead of competing for the Fermi surface), the su- 
perfiuid and the CDW coexist. Using Green's functions, 
we derive the exact critical repulsion Vc in two limiting 
cases; we find Vc = 4^2 for the two-particle system and 
Vc = 2^/21-2 for the half-filled system. Using finite size 
scaling analysis [with Vc{N) — Vc{oo) ex 1/A^^ for system 
size N] , we also obtain the critical repulsion at other fill- 
ings numerically by using modified Lanczos algorithmic. 

Our t2 — V model can be mapped onto an 
extremely anisotropic Heisenberg model (with 
next-nearest-neighbor XY interaction and nearest- 



neighbor ising interaction) and is of the form 
J^ E.iStS-+2 + H.c.) + J,| E, SfSf^,. While Heisen- 
berg model was amenable to solution through the Bethe 
ansatz, the addition of next- nearest- neighbor interaction 
requires an alternate route for its solution. Thus, the 
Majumdar-Ghosh modeli^ {JT.^[S^■Si+l + \Si-S^+2\), 
which is a special case of the isotropic Heisenberg model 
with nearest- and next-nearest-neighbor coupling, can 
be solved by an alternate elegant approach yielding 
an exact ground state cast in a simple valence-bond 
form. Our spin model too lends itself to exact instability 
solutions (by the Green's function method) in the two 
limiting cases of two-magnons and antifcrromagnetic 
ground state. 
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II. NUMERICAL STUDY OF THE t2-V MODEL 



We begin by identifying the Hamiltonian of the t2 
model for hard-core-bosons (HCB). 



V 



N 



TV 
+^"^^1^71^ + 1, (1) 

where bj is destruction operator for a HCB, V > 0, 
Uj ~ bjbj, and N (which is even) is the total number of 
sites. We will study numerically the quantum phase tran- 
sition (QPT) in the t2 — V model by assuming periodic 
boundary conditions. We will characterize the transition 
through the structure factor and the superfluid density. 
The structure factor is given by 



N 



5(fc) = ^e^*^'W^(Z), 



(2) 



(=1 



where W{1) is the two-point correlation function for den- 
sity fluctuations of HCBs at a distance I apart (when 
lattice constant is set to unity). 



N 



W{l)^^J2i{n,n,+i)-{n,){n,+i)] 



(3) 



f 



wavevector k = 2^ with n=l,2, ,N; filling-fraction 

rij) ~ -jf with Np being the total number of HCBs 
in the system. Now, we observe that N,,, ~ ^ ■ Uj 
(No = J2j ^j) ^^^ number operator which gives the 
total number of electrons at even (odd) sites commutes 
with the Hamiltonian; hence the following simple expres- 
sion resultsi 



Si-K) 



^{{N, - Nof) A{N, - Nof 



N 



N 



(4) 



FIG. 1. (Color online) Plots of rescaled structure factor S*(7r) 
and superfluid fraction Us at various filling factors f obtained 
using modified Lanczos technique. The calculations were per- 
formed at fiUings f = 1/2, 1/4 with system size A'^ = 16 while 
at f = 1/3 size A*' = 12 was used. At a critical repulsion there 
is a striking discontinuous transition in both ^'(Tr) and n^; 
while S* (tt) jumps from its minimum to maximum, there is a 
significant drop in n^. 



where N^ {No) are the number of HCBs in even (odd) 
sites. Thus the minimum value of ^(Tr) = corresponds 
to equal number of particles in both the sub-lattices 
whereas the maximum is given by 



[S{^)] 



ANl 

N 



(5) 



indicating a single sublattice occupancy. To study the 
QPT, we rescale the value of S'(7r) as follows: 

S{^) 



S*{^) 



[S[^)]r 



(6) 



with S* (tt) representing the order parameter that varies 
from to 1 during the phase transition. 

Next, we will outline our procedure for calculating the 
superfluid density by threading the chain with an in- 
finitesimal magnetic flux 9. The superfluid fraction is 
given by22i2i 

7V2 \ld'^E{e)' 
2" 



Nr,t 



p^<^Sf 



00^ 



(7) 



9=0 



where E{9) is the total energy when threaded by flux 
6 and t^s = H'^/2m is the efi^ective hopping term which 
for our t2 — V model is given by toff = 4^2- The total 
energy for the case V = when threaded by a flux 6, is 
expressed as 



E{0) 



-2i2^cos[2(fc 

k 



e/N)] 



(8) 



When, both sub-lattices are occupied, the superfluid frac- 
tion is given by2^ 



(^) 



Nr, 



\ N , 



(9) 



whereas when only one sub-lattice is occupied it is ex- 
pressed as 



1 



sm 






N„ 



-n m 



(10) 



where periodic boundary conditions have been taken. 

Here, we do not calculate the Bose-Einstein conden- 
sate (BEC) occupation number ng because, for a sys- 
tem of HCB in a one-dimensional tight-binding lattice, it 
varies as C{l)\N in the thermodynamic limit with the 
coefficient C(f) depending on filling fi-action f— i2i; con- 
sequently, the condensate fraction rio/Np ex 1/vN — >■ 0. 
Next, in the presence of repulsion (as argued in Ref. l22l ). 
we expect the BEC occupation number ng to again scale 
as VN', however, the coefficient of y/N will be smaller 
due to the restriction on hopping imposed by repulsion. 

When rings with even number of sites were used, at 
all fillings f, we found that the order parameter S*{tt) 
jumps from to 1 at a critical value of repulsion Vc in- 
dicating that the system transits from equally populated 
sub- lattices (i.e., ising symmetry) case to a single sub- 
lattice occupancy, i.e., a period-doubling CDW state [see 
Fig. ([1])]. Concomitantly, as can be seen from Fig. ([1]), 
there is a sudden drop in the superfluid fraction n^ at the 
same critical repulsion. At half-filling, where the super- 
fiuid fraction vanishes above a critical repulsion because a 
single sub-lattice is completely filled, the transition shows 
that superfluidity and CDW state are mutually exclu- 
sive. On the other hand, at all non-half-flllings, we see 
that the system undergoes a QPT from a superfluid to 
a supersolid (i.e., a homogeneously coexisting superfluid 
and CDW) state. In Fig. ([T]), it is of interest to note 
that the values of n^ at y = and V > Vc are exactly 
those predicted by Eq. ([9|) and Eq. pO|) respectively. 

Next, using the finite size scaling expression (as de- 
rived in Appendix A) relating Vc{N) (critical repulsion 
at a finite N) to Vc{oo) (critical repulsion of an infinite 
system) 



v.iN) - Vci^) = ^ + |r 



(11) 



we obtain the critical Vc(oo) depicted in Fig. ([2]) (from 
finite size systems calculations) at various fillings. At 
half-filling, from finite size scaling analysis we obtain 
that Vc{oo) sa 2.83; in the next section, we show (using 
Green's function analysis) the exact result Vc{oo) = 2-\/2. 
For systems with 2 HCB and iV = 4, 6, 8, 10, 12, 14, 16, 18, 
and 20, we find numerically that Vc ~ 4.00; in the next 
section, we obtain the exact result (using Green's func- 
tions) that Vc{N) = 4 for any system size A^ > 4. We 
will now cast the t2 — V model (for HCBs) as an ex- 
tremely anisotropic Heisenberg spin model by identifying 
S^ = b\ S^ = b, and S^ = n — ^. The resulting spin 
Hamiltonian is of the form 



N 



t2^istA 



i+1 



H.c. 



N 

vY^s^s^,. 




FIG. 2. (Color online) Plot of Vc{oo) (critical repulsion for an 
infinite system) obtained from Green's function analysis for 
half-filled (f = 1/2) and two HCB systems (f — >■ 0) and using 
finite size scaling at various other fillings f . 



We see that the above spin model has NNN interaction in 
the transverse (XY) channel while it has NN interaction 
in the longitudinal (Z) channel. The energies at various 
fillings Np/N for the t2~V model correspond to various 
magnetizations {N — 2Np)/N for the spin model. From 
a plot of the energies at various magnetizations of the 
spin model, as depicted in Fig. ^, we see that the en- 
ergy increases with magnetization with the ground state 
corresponding to zero magnetization. From the fact that 
always the critical repulsion Vc < 4, it should be clear 
that the energy at all fillings and system sizes for Vc > 4: 
is obtained from a tight-binding model with Np particles 
in one sub-lattice only. Thus at Vc > 4, the energy will 
certainly increase with magnetization. 




(12) 



FIG. 3. (Color online) Plots of energy Eg (obtained using 
modified Lanczos in a system with A'^ = 16 sites) for the ex- 
tremely anisotropic NNN Heisenberg model at various mag- 
netizations {N — 2Np)/N corresponding to fillings f = Np/N 
for the ti — V model. 



III. ANALYSIS OF THE HALF-FILLED t2 - V 
MODEL 

We study the following t2 — V Hamiltonian in rings 
with even number of sites by considering two sublattices 
C and D and using periodic boundary conditions: 

N/2 

ff=-t2^(cJc,+i+H.c.) 

i=l 
N/2 

-t2^(4'^»+i+H.c.) 

i=l 
N/2 

+F^4d,(cjc. + cj_ic,-i), (13) 

i=l 

where c and d denote destruction operators of HCBs 
in sublattices C and D, respectively. To understand 
the striking discontinuous phase transition at half-filling, 
(i.e., the transition from equal occupation of both sublat- 
tices to occupation of only one sublattice at a critical Vc) 
we begin by recognizing that when the system is on the 
verge of completing the phase transition, the system will 
pass through the state where there is one particle in one 
sublattice and one hole in the other sublattice. Hence, 
we now consider instability for the case of one particle 
in sublattice C and one hole [with destruction operator 
denoted by h (= d^)] in sublattice D; we then rearrange 
the above equation as 

N/2 



H=-t2j2{4c,+i+R.c.) 

i=l 
N/2 

+t2j2{hlh,+i+'ti.c.) 

1=1 

N/2 
-VY^{h\K ~ l){c\c^ + cUc,-l). (14) 



We define the particle-hole Green's function as foUowa^^: 

g^' = ^(fc,0|G(w)|fc,n)^, (15) 

where 

1 



GH 



[uj + iri-HY 
and the particle-hole state \k,n)h 

N/2 



(16) 



N/2 

k,n),^—Y,e^''i'''+^)4hljO), (17) 



with k being the total momentum of the particle-hole 
system, na the separation between the particle and hole, 
and a the lattice constant of a sublattice. Using the con- 
dition 



for n ~ 0, we obtain 

(w + iv)go - 1 + h{k, 0\G{uj){H)\k,0)i, 
= 1 + Vg^ - m-i €mikal2)g\ 

+i2t2Sin{kal2)g\. (19) 

From Eq. (fT8|) . for n = 1. we get 

{uj + ii^)g'l = h{kMG{u^){H)\k,l)h 
= Vg'l - i2t2 sm{ka/2)g^ 

+i2t2 sm{ka/2)g^. (20) 

Similarly, for ?i 7^ 0, 1, we derive 

{cj + i,j)gl[^^,{k,0\G{oj){H)\k,n)h 

= 2^5^ - i2t2 sin(fca/2).g,'j+i 

-Fi2t2sin(fca/2)5r^_i. (21) 

As V increases to the critical interaction Vc, the energy 
given by Eq. ([T3l) becomes (i.e., energy of all HCBs 
filling up one sublattice). This would correspond to the 
level-crossing instability [in the system represented by 
Eq. (fT4|) ] where the particle HCB quits its sublattice 
and goes into the sublattice of the hole. 

Now, let 



uj + iT]-2V _ 1 



z, 



(22) 



where Fk = i2t2Sm{ka/2). Then, Eq. (gT]) takes the 
simple form 



1 



Z ] 9n ~ 9n~l 5n+li 



whose solution is of the form 



9n 



atz"-' 



^t 



{-zY 



(23) 



(24) 



where af {o-i) and /3^ (/3f ) correspond to n > 1 (n < 0). 
The transition occurs at the critical value of V that makes 
the overall energy 0. Let V = 2t27sin(fca/2). The overall 
energy is less than ~it2+2V (for V > 0). It is important 
to note that, for the groundstate, k = n/a for any V. 
This can be seen by first noting that when V = 0, total 
momentum in the minimum energy state is 7r/a; next, 
turning on V does not change the total momentum. Then 
for k = ir/a, to get overall energy to be 0, we need the 
inequality 7 > 1. To get the instability condition, we 
consider the case w = so that the Green's function ^q 
diverges at the transition point. It then follows from Eq. 
(HH) that 



i27= - 
z 



Zi 



(25) 



which implies that 



<5o,„ = h{k, G\G{lo){u: + irj- H)\k, n)h, (18) 



z = J(-7+ Vt-1), 



(26) 



and hence \z\ < 1 for 7 > 1. Let us first consider tlie 
case n > 1. From tlie above Eq. (|24p . it is clear that, for 
n -^- 00. 5^ is finite only for /3^ ~ 0. Thus for n > 1, 



9n+l — ^9m 



(27) 



which implies that 53 = zg!^; then, from Eq. (|23p it 
follows that 



9^ = ^5^ 



(28) 

Next, for the case n < 0, we see that g^ is finite, for 
n — > —00, only when aj~ = 0. Thus for n < 0, 



5,ti 



-Z9n, 



(29) 

which implies that g'ij = —zg^^] then from Eq. (P^ we 
obtain 



5-1 



'^ffo- 



From Eqs. dH]), d^H]), dUl), and (gni), we obtain 

1 



50^ 



(w + ir; - y + zFfc) + (^+,^„{>+,j.^) 



and 



.9i 



.90'^^ 



[uj + i-q-V + zFk) ' 



(30) 



(31) 



(32) 



where uj = Q, V = 2^27, z = i(— 7 + ^7^ — 1), and 
Ffc = 12^2- It then follows that ^q diverges when {V — 
zFkY + Fl = 0, i.e., when 7 = \/2. Thus the instability 
condition is Vc = 2-\/2i2- (We now see that the total 
energy at transition is indeed less than —4^2 + 2V). It is 
important to note [as can be seen from Eq. p2p ] that, 
when 90 diverges, g^ also diverges; consequently, from 
Eqs. dill), dSni), and dSO]) we see that all 5^ diverge (i.e., 
even when n > 1 and n < 0). 



IV. TWO HCBS llSi t2-V MODEL 

We study the t2 — V Hamiltonian [described by Eq. 
(|13p ] in rings with even number of sites and for periodic 
boundary conditions. We consider one HCB in sublattice 
C and one HCB in sublattice D; the corresponding two- 
particle Green's function is defined by 



.g„=(fc,0|G(c^)|fc,n), 



(33) 



where G(w) is given by Eq. (J16p and the two-particle 
state |fc,n) is defined as 

N/2 

%n) = ^e'K«.+¥)cJdt^jO), (34) 



with k representing the total momentum of the two- 
particle system. Then, the following equations hold for 
the Green's functions g„: 

(lu + irj — T^),go = 1 ^ 2^2 cos{ka/2)gi 

-2f2Cos(fca/2)g_i, (35) 



{uj + ir] - V)gi 

and for n 7^ 0, 1 

(a; -I- iri)gn = - 



-2^2 cos{ka/2)g2 
-2t2cos{ka/2)go, (36) 



-2^2 cos(fca/2)5„+i 

-2^2 cos(fca/2)g„_i. (37) 



As V increases to the critical Vc, the energy given by 
Eq. dEI becomes — 4*2 cos(27r/A^) (i.e., the minimum 
energy of 2 HCBs in the same sublattice); this would 
correspond to the instability for Eq. dI31) where one HCB 
quits its sublattice and goes into the sublattice of the 
other particle. Here, we make the key observation that 
fc = for the groundstate at any V; to understand this, 
we first note for T^ = 0, the total momentum is zero in the 
minimum energy state; next, we recognize that turning 
on V does not change the total momentum. 

Now, to obtain the instability , we take 



w _ 4*2 cos(27r/A^) 
(-2*2) ~ 2^2 



2cos(27r/7V) 



z + 



(38) 



We set z = e^'^'^l^ and assume r^ — > faster than 
sin(27r/A^) as iV — ?► 00 in the thermodynamic limit. We 
also take F/(2i2) = 27. Then Eqs. dMI, dMl), and dSZ]) 
become 



[(z + 1/z) + 27]5o = -1/(2*2) +51+5- 



(39) 



(40) 



[(z + 1/z) + 27].gi = .92 + .90, 

and for n 7^ 0, 1 

[(z + l/z)]g„ = 5„+i + 5„_i. (41) 

Without loss of generality, we assume 

51 =a2Z-F/32/z, (42) 

and 

g2=a2z2 + /32/z'. (43) 

Then using Eqs. dH]), d32)' and dSl), we obtain for n = 
3,4,..., iV/2 the expression 



g„ = q;2z" -f /32/z' 



(44) 



It is important to note that [zl = 1 and hence the Green's 
functions do not decay with n\ From the fact that at two- 
particle momentum fc = 0, |/s, n) = |/s, n — -/V/2), we see 
that gjv = (70 and gN__i — g-i- Then, using Eq. (|44p 

and the relation z = e'^'^/^, we get 

5-1 = 5f-l = "2Z 2 +P2lZ''- 

= -{a2lz + hz), (45) 



and 
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= -(a2 4-/32). 



(46) 



We are now ready to solve for the Green's functions g„ 
using Eqs. dSHD, &, ©, dSD, (US)), and dH]). We get 
the foUowing equations: 



a2[z + 7] + P2[l/z + 7] - 1/(4*2), 



(47) 



and 



a2[l + jz]+P2[l+l/z] = 0. (48) 

It then foUows from the above two equations that 

z + 7 



a2 



4*2(^2 +72 - 1 -72^2)^ 



(49) 



which diverges when 7 = ±1. We choose 7 = 1 for 
repulsive V. Furthermore, 



/?2 = -"2 



1 +7Z 
1+7/2' 



(50) 



which for 7 = 1 yields /32 = —0:22. Thus we see that 
go = —{a2 + P2) = — 012(1 — 2) diverges for 7 = 1 or 
Vc = 4*2 ■ We also find that for 7^ n < 



N 



5n = "22 +P2/2 = ^2 I e ' -e ^ ^' 1,(51) 

also diverges since 2n— 1 f^ A'^. T/ie instability/divergence 
condition Vc = 4*2 *s independent of N and hence is valid 
in the thermodynamic limit as well! Another interesting 
observation based on /32 = —022 is that 



g-k = .94-fc = -012/ z^ - I32Z^ 

= /32/2'=+l + 022* 

= 5fc+i- 
V. CONCLUSIONS 



(52) 



In the t2 — V model, for V » t2, energy consider- 
ations show that the system always has only one sub- 
lattice populated for less than half-fillingi. Addition- 
ally, we demonstrate that our t2 ~ V model undergoes 
a dramatic first-order QPT at non-half-fillings from a su- 
perfluid to a supersolid state with homogeneously coex- 
isting U(l) symmetry broken superfluid and a period- 
doubling CDW. Using two-particle Green's functions we 
derive the QPT condition exactly for the half-filled case 
and the instability condition for the 2 HCB case for all 
system sizes. The model studied is a limiting case of the 
strong electron-phonon coupling in one-dimensional sys- 
tems with cooperative breathing mode. A similar analy- 
sis in two-dimensions is left for future studies. 
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Appendix A: Finite size scaling analysis for HCBs in 
t2 — V model 

In this appendix, we will outline our approach to car- 
rying out finite size scaling analysis. Let us first consider 
a non-interacting system (with A^ sites and A'j, particles) 
described by a tight-binding Hamiltonian. For even num- 
ber of particles, the ground state has particles occupying 
momenta "]v^ with —Np/2 < m < Np/2— 1; whereas 
for odd number of particles the ground state is repre- 
sented by particle momenta ^ ™j^ with —{Np — l)/2 < 
m < {Np — l)/2. Thus the ground state wavefunction 
(due to the occupied momenta) is an even function of 
1/A^. Hence, for the case corresponding to after the phase 
transition where all the particles are in the same sub- 
lattice C, the energy of the ground state \<j>o) is given 
by 

N/2 

Ei = iZ^M- t2{cU+i + H.c.)|0o>, (Al) 



where c is the destruction operator for a HCB in sub- 
lattice C. Upon taking into account discrete translational 
symmetry, we get 



2Ei 



/)o|-t2(clc,+i+H.c.)|0o)- (A2) 



Since |0o) is even in 1/iV, we note that -^ is also even 
in l/N. 

Next, consider the interacting system characterized by 
the following t2~V Hamiltonian in rings with even num- 
ber of sites (iV) and with periodic boundary conditions: 

N/2 

i?=-*2^(clc,+i+H.c.) 
j=l 

N/2 

-t2^{dld,+i+il.c.) 

i=l 

N/2 

+V^dld^iclc^ + ct_iQ-i), (A3) 

i=l 

where c {d) is the destruction operator for HCB in sub- 
lattice C (I?)and V > 0. Upon invoking refiection sym- 
metry, we note that the ground state will be invariant 
when the sign of momenta is reversed; equivalently the 
ground state \iPq) is an even function of 1/A'^. The ground 



state energy, before the phase transition, is given by 

N/2 

Eii = Y.^M - i2[(clc,+i +H.C.) + (4d,+i +H.c.)]|V'o> 

i=l 

N/2 

+ Y,{MVd\d,{c\c, + 4_iC,„i)|Vo)- (A4) 

Upon recognizing discrete translational invariance, we see 
that 



2Eu 

N 



= (^ol - t2[{c\c,+i + H.c.) + {d\d,+i + H.c.)]|V'o) 



-i^^.(J..^J 



{'^pQ\Vd]d^{clci + c!_iCj_i)|-0o) 



(A5) 



Since, \'iPq) is even in 1/A^, it follows that -^ is also even 

in 1/iV. 



Now, at the transition point (corresponding to a criti- 
cal interaction 14), -^ — -^ = 0; this implies that 

(V'o|i2[(cIc,+i +H.C.) + (d|d,+i +H.c.)]|^o) 



K 



{i})o\d\di{c\ci + c\_-^Ci-i)\ilJa) 
(0o|t2(4c,+i+H.c.)|0o) 
{iiQ\d\di{c\ci + c\_^Ci-i)\iiQ) 



(A6) 



In the above equation, because both numerator and de- 
nominator of all the terms on the right-hand side are even 
in 1/-/V, it follows that Vc is also even in \/N . 

It is also important to note that we used general argu- 
ments to show that the ground state energy ^ff- is even in 
1/A'^; these arguments can be extended to show that the 
ground state energies of other interacting systems such 
as the t — V model are also even functions of 1/iV. 
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